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HIGHLIGHTS 


►  Analytic  functions  describing  the  battery  model  parameters  are  optimized. 

►  EKF  based  on  the  optimal  analytic  model  is  adopted  as  the  SOC  estimator. 

►  The  robustness  of  the  SOC  estimator  against  varying  loading  profiles  is  evaluated. 

►  The  robustness  of  the  SOC  estimator  against  varying  temperatures  is  analyzed. 

►  The  robustness  of  the  SOC  estimator  against  varying  aging  levels  is  assessed. 
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Battery  State  of  Charge  (SOC)  estimation  is  an  important  function  for  battery  management  systems  and 
critical  for  the  reliable  operations  of  batteries.  This  paper  analyzes  the  robustness  of  SOC  estimation 
algorithms  for  two  types  of  Li-ion  batteries  under  varying  loading  conditions,  temperatures  and  aging 
levels.  Based  on  the  model  templates  identified  in  an  earlier  research,  the  model  parameters  are 
determined.  The  Extended  Kalman  Filter  (EKF)  technique  is  then  adopted  as  the  SOC  estimation  algo¬ 
rithm.  The  robustness  of  the  estimator  against  varying  loading  profiles  and  temperatures  is  evaluated 
and  compared  against  the  Coulomb  counting  method.  We  subsequently  used  data  from  cells  that  have 
significantly  aged  to  assess  the  robustness  of  the  SOC  estimation  algorithm.  Finally,  the  need  for  model 
parameter  updates  is  analyzed. 

©  2012  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Electrified  vehicles  including  battery  electric  vehicles  (BEVs), 
hybrid  electric  vehicles  (HEVs)  and  plug-in  hybrid  electric  vehicles 
(PHEVs)  have  been  actively  studied  and  developed.  A  critical 
element  for  their  successful  commercialization  is  technologies  for 
reliable  battery  operations.  Battery  management  systems  (BMS) 
have  been  designed  to  provide  monitoring,  diagnosis,  and  control 
functions  to  enhance  the  operations  of  battery  packs.  A  key  func¬ 
tion  of  BMS  is  to  accurately  estimate  battery  state  of  charge  (SOC). 
Poor  SOC  estimation  can  result  in  larger  SOC  swing  than  specified, 
and  can  lead  to  reduced  cycle  life  or  lower  efficiency. 
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The  Coulomb  counting  method  was  often  used  as  a  core  tech¬ 
nology  for  battery  SOC  estimation  [1—3].  This  method  is  easy  to 
implement  but  has  three  challenges.  First,  the  initial  SOC  at  key-on 
must  be  estimated  accurately  because  subsequent  estimates  may 
be  biased  by  the  initial  SOC  error  [4].  Secondly,  the  method  is  highly 
dependent  on  the  accuracy  of  the  current  sensor.  The  current 
sensor  is  often  subject  to  noise,  drift  induced  by  temperature  and 
other  uncertainties.  Finally,  the  battery  capacity  reduces  with  age 
and  will  affect  SOC  calculation.  To  remedy  these  three  challenges, 
periodic  resets  are  needed.  In  actual  vehicle  operations,  the  open 
circuit  voltage  (OCV)  is  often  used  to  calculate  initial  SOC.  Since 
there  is  typically  a  monotonic  relationship  between  OCV  and  SOC. 
Inverting  this  algebraic  relation  leads  to  a  SOC  estimate.  However, 
error  can  exist  if  the  battery  is  not  fully  relaxed.  Sometimes  the 
OCV— SOC  curve  is  very  flat  in  the  middle  region  (e.g.,  between  30 
and  80%  SOC),  particularly  for  lithium  iron  phosphate  cells.  A  small 
OCV  error  can  result  in  a  large  SOC  error.  Hence,  the  OCV-based 


0378-7753 /$  -  see  front  matter  ©  2012  Elsevier  B.V.  All  rights  reserved. 
http://dx.doi.Org/10.1016/j.jpowsour.2012.06.005 


210 


X.  Hu  et  al.  /  Journal  of  Power  Sources  217  (2012)  209-219 


Table  1 

Equations  of  the  best  performing  models  identified  in  Ref.  [24]. 

Model  Equations 

l4+i  =  exp(— At/Ti)t4  +  Ri[l  -exp(-At/T!)]/k 
Vk  =  VOc(zk)-R0Ik-Uk 

where  I,  V  and  z  are  current,  output  voltage 
and  SOC,  respectively.  At  is  the  sampling  time 
and  V0c  is  OCV.  R0  is  the  internal  ohmic 
resistance  which  depends  on  the  current 
direction.  U  and  ii  =  RiCi  are  the  voltage  and 
time  constant  of  the  RC  network. 

/  Uk+1  =  exp(-At/i1)t/k+R1[l  -exp(-At/T-,)]4 
\  hk+ 1  =  exp(-|/c4At|)hk  +  [1  -  exp(-|/c/kAt|)]H 
Vk  =  Voc(Zk)  -  Rohi  -  Uk  +  hk 
where  h  is  the  hysteresis  voltage,  k  is  a  decaying 
factor  and  H  is  the  maximum  amount  of 
hysteresis  voltage  which  is  positive  for  charge 
and  negative  for  discharge. 


reset  is  not  suitable  for  batteries  that  are  mostly  half-charged.  For 
EVs/PHEVs,  the  reset  can  be  more  effectively  done,  since  the  SOC 
frequently  reach  the  two  ends  where  the  OCV-SOC  slope  is  steep. 

Many  SOC  estimation  models  have  been  proposed,  such  as 
artificial  neural  networks  based  models  [5—7],  fuzzy  logic  models 
[8,9]  and  support  vector  regression  (SVR)  based  models  [10,11  ].  The 
robustness  of  these  models  strongly  relies  on  the  quantity  and 
quality  of  the  training  data  set.  A  limited  training  data  set  (e.g., 
obtained  from  new  cells)  may  result  in  limited  model  robustness, 
thus  reducing  the  applicability  of  the  model. 

Kalman  filter  and  other  observer-based  approaches  have  also 
been  used  to  estimate  the  battery  SOC.  These  methods  use  output 
feedback  and  can  have  better  robustness  than  non-feedback 
methods.  In  Refs.  [12-14],  an  Extended  Kalman  filter  (EKF),  based 
on  several  nonlinear  state-space  models,  was  used  to  estimate  the 
SOC  of  an  HEV  Li-polymer  cell.  A  central  difference  Kalman  filter 
(CDKF)  using  a  nonlinear  enhanced  self-correcting  battery  model 
was  also  developed  to  estimate  SOC  [15,16].  An  unscented  Kalman 
filter  (UKF)  by  means  of  a  nonlinear  electrochemical  battery  model 
was  also  used  to  estimate  SOC  of  a  lithium-ion  cell  [17].  The 
performance  of  these  filters  depends  on  the  model  accuracy.  The 
sigma-point  based  CDKF/UKF  often  provides  better  estimates  than 
EKF  at  the  expense  of  higher  complexity  and  computational  cost. 
Several  variants  of  Kalman  filter  were  also  studied  with  similar 
results  [18-20].  Besides  Kalman  filters,  sliding-mode  observers 
[21,22]  and  linear  parameter- varying  observers  [23]  were  also 
applied  to  predict  the  battery  SOC. 

Most  of  the  estimation  methods  described  above  were  validated 
using  battery  data  under  a  narrow  set  of  scenarios,  without 


Table  2 

Candidate  OCV  functions  evaluated  in  this  paper. 


No. 

OCV  functions  and  description 

(1) 

Voc  =  /Co  -  /q/z  -  I<2z  +  I<3  In  (z)  +  JC4  In  (1  -  z)  [13] 
where  VDc  and  z  are  the  OCV  and  SOC,  respectively. 

The  optimization  variable  vector  6  =  [/C0,  K\,  /C2,  /C3,  JC4]. 

(2) 

Voc  =  K0  +  aAz  +  a2(\  -e~^z)  +  a4(\  -  e-a^~z)  [26] 
where  6  =  [JC0,  a3,  a2,  a3,  a4,  a5]. 

(3) 

Vac  =  Ko-cti/z+^e-^'-*)  [27] 
where  0  =  [ I<0 ,  a3,  a2,  a3]. 

(4) 

Voc  =  a^e~a2Z  +  a3  +  a4z  +  a5z 2  +  a6z3  [28] 
where  6  =  [ai,  a2,  a3,  a4,  as,  a^]. 

(5) 

Voc  =  ol\Z6  +  a2z5  +  a3z4  +  a^3  +  a&2  +  a^z  +  a7  [29] 
where  6  =  [a\,  a2,  a3,  a4,  as,  ai6,  a2]. 

SOC 

Fig.  2.  OCV  fitting  results  for  LiNMC  cell  at  22  °C. 

exploring  different  temperatures,  battery  ages,  or  highly  transient 
loadings.  In  other  words,  the  robustness  of  these  SOC  estimation 
algorithms  was  not  sufficiently  assessed.  For  example,  many  SOC 
estimation  approaches  mentioned  above  were  evaluated  under 
only  one  battery  loading  profile  and  one  environmental  tempera¬ 
ture.  Moreover,  the  performance  and  robustness  of  these  SOC 
algorithms  against  aging  were  not  adequately  studied.  A  key 
contribution  of  this  paper  is  that  the  SOC  estimator  performance  for 
two  types  of  Li-ion  batteries  was  evaluated  under  different  loading 
profiles,  temperatures,  and  cell  aging  levels. 
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Fig.  1.  The  test  schedule  to  collect  battery  cell  data. 


Fig.  3.  OCV  fitting  results  for  LiFeP04  cell  at  22  °C. 
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Three  candidate  analytic  functions  for  RC  and  hysteresis  parameters. 
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Type 

Functions  and  description 

(1)  Polynomial  [29] 

/=CnZn  +  Cn_iZn-1+"-+C0 

where  /  and  z  are  model  parameter  and  SOC, 
respectively,  n  is  the  degree  of  the  polynomial. 
cn  is  a  coefficient  to  be  optimized. 

(2)  Exponential  [28] 

/  =  c\  e~C2Z  +  c3 

(3)  Power  series  [30] 

/  =  Ci  zC2  +  c3 

Based  on  the  preferred  model  structures  identified  in  our  earlier 
work  [24],  the  optimal  analytic  functions  explicitly  depicting  the 
dependency  of  model  parameters  on  SOC  and  temperature  are 
determined.  These  models  are  then  used  in  EKF  to  estimate  SOC. 
The  robustness  of  the  estimator  under  varying  loading  profiles, 
temperatures  and  cell  age  are  then  analyzed. 

The  remainder  of  this  paper  is  organized  as  follows.  In  Section  2, 
determination  of  the  optimal  analytic  functions  to  explicitly 
delineate  the  model  dependence  on  SOC  and  temperature  is 
described.  In  Section  3,  the  EKF-based  SOC  estimator  is  introduced. 
The  evaluation  results  of  the  estimator  robustness  are  discussed  in 
Section  4,  followed  by  conclusions  in  Section  5. 

2.  Optimization  of  model  dependency  on  SOC  and 
temperature 

In  Ref.  [24],  we  studied  twelve  equivalent  circuit  battery 
model  structures  reported  in  the  literature,  and  compared  their 
complexity,  training  and  validation  accuracies  using  test  data 
from  multiple  cells.  The  preferred  model  structures  for  the 
lithium  nickel-manganese-cobalt  oxide  (LiNMC)  UR14650P  and 
lithium  iron  phosphate  (LiFeP04)  APR18650M1A  cells  were  found 
to  be  the  first-order  resistance-capacitance  (RC)  model  and  the 
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Fig.  4.  Schematic  diagram  of  the  two-step  optimization  process. 


C  FUDS  data 

Fig.  5.  Comparison  results  of  polynomials  with  different  degrees  for  the  LiNMC  cell. 


first-order  RC  model  with  one-state  hysteresis,  respectively.  The 
equations  of  the  two  models  are  shown  in  Table  1.  The  aim  of  Ref. 
[24]  was  to  find  the  optimal  model  structures  and  the  parameters 
of  all  the  twelve  models  were  assumed  to  be  constant  (inde¬ 
pendent  of  SOC).  To  explore  the  full  potential  of  these  two 
preferred  model  structures,  in  this  paper  we  will  establish  the 
dependence  of  model  parameters  on  the  battery  SOC  and 
temperature.  The  test  schedules  shown  in  Fig.  1  are  designed  to 
excite  and  age  the  two  types  of  battery  cells.  The  datasets  are 
described  in  details  in  Ref.  [24].  The  three  dynamic  test  cycles 
(Hybrid  Pulse  Test  DST  and  FUDS)  are  used  to  identify  the  SOC 
and  temperature-dependent  analytic  functions. 
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Table  4 

Optimal  coefficients  of  the  analytic  functions  for  the  LiNMC  cell. 


Parameter  Coefficients 


OCV 


Charging  R ^ 


Discharging  Rj 


R i  (RC  network) 


Time  constant  ii 
(RC  network) 


Voc  =  a\ z6  +  a2z5  +  a3z4  +  a. 4Z3  +  a5z2  +  a^z  +  a7  (0.1  <  z  <  0.9) 
where 

=  {  -0.1674(Te  -  10)  +  2.9192,  10  <  Te  <  22  (Te  :  Temperature);  0.1939(Te  -  22)  +  0.9109,  22  <  Te  <  35} 
a2  =  {  0.3394(Te  -  10)  -  5.8608,  10  <  Te  <  22;  -0.3675(Te  -  22)  -  1.7875,  22  <  Te  <  35} 

a3  =  {  -0.1811(Te  -  10)  +  1.3408,  10  <  Te  <  22;  0.0469(Te  -  22)  -  0.8330,  22  <  Te  <  35} 

aA  =  {  -0.0614(Te  -  10)  +  5.1684,  10  <  Te  <  22;  0.3023(Te  -  22)  +  4.4320,  22  <  Te  <  35} 

a5  =  {  0.0948(Te  -  10)  -  4.4983,  10  <  Te  <  22;  -0.2336(Te  -  22)  -  3.3606,  22  <  Te  <  35} 

a6  =  {  -0.0290(Te  -  10)  +  1.5815, 10  <  Te  <  22;  0.0648(Te  -  22)  +  1.2331,  22  <  Te  <  35} 
a7  =  {  0.0025(Te  -  10)  +  3.5600,  10  <  Te  <  22;  -0.0063(Te  -  22)  +  3.5901 ,  22  <  Te  <  35} 

Rq  =  cxe~C2Z  +  c3  (0.1  <  z  <  0.9) 
where 

Cl  -  { -0.01 16(Te-  10) +0.1728,  10  <  Te  <  22;  0.0741  (Te  -  22)  +  0.0337,  22  <  Te  <  35} 
c2  =  {  0.5632(Te  -  10)  -  24.5426,  10  <  Te  <  22;  -1.0925(Te  -  22)  -  17.7848,  22  <  Te  <  35} 
c3  =  {  -0.0009(Te  -  10)  +  0.1086, 10  <  Te  <  22;  -0.0005(Te  -  22)  +  0.0972,  22  <  Te  <  35} 

R+  =  ci  e~C2Z  +  c3  (0.1  <  z  <  0.9) 
where 

C\  =  (0.0072(Te- 10) +  0.0761,  10  <  Te  <  22;  -0.0125(Te  -  22)  +  0.1621,  22  <  Te  <  35} 
c2  -  {  -2.0975(Te  -  10)  -  4.7917,  10  <  Te  <  22;  -0.6137(Te  -  22)  -  29.9621,  22  <  Te  <  35} 
c3  =  {-0.0010(Te- 10) +0.1142,  10  <  Te  <  22;  -0.0006(Te  -  22)  +  0.1021,  22  <  Te  <  35} 

R\  =  Cie_C2Z  +  c3  (0.1  <  z  <  0.9) 
where 

Ci  =  (0,  10  <  Te  <  22;  0,  22  <  Te  <  35}  c2  =  {  -0.1403(Te  -  10)  -  21.9590,  10  <  Te  <  22;  -0.6098(Te  -  22)  -  23.6430,  22  <  Te  <  35} 
c3  =  {-0.0015(Te-  10) +0.0608,  10  <  Te  <  22;  -0.001 1  (Te  -  22)  +  0.0433,  22  <  Te  <  35} 
ii  =  Cie_C2Z  +  c3  (0.1  <  z  <  0.9) 
where 

Ci  -  (2.1871(Te- 10) +  19.3967,  10  <  Te  <  22;  -1.9460(Te  -  22)  +  45.6424,  22  <  Te  <  35} 
c2  =  {  -4.1106(Te  -  10)  -  0.6723,  10  <  Te  <  22;  1.9550(Te  -  22)  -  50.0000,  22  <  Te  <  35} 
c3  -  {0.5091(Te- 10) +  32.7608,  10  <  Te  <  22;  1. 7591  (Te  -  22)  +  38.8700,  22  <  Te  <  35} 


2.1.  OCV  function  identification 

Compared  to  lookup  tables,  analytic  OCV-SOC  functions  can  be 
beneficial.  For  example,  function  derivatives  can  be  calculated  more 
easily  and  accurately.  Several  candidate  functions  proposed  in  the 
literature  to  depict  the  battery  OCV  are  summarized  in  Table  2.  We 
fit  the  five  candidate  functions  shown  in  Table  2  to  the  12  OCV 
points  obtained  in  [24]  by  nonlinear  least-squares  optimization. 
The  results  of  the  LiNMC  and  LiFeP04  cells  at  temperature  =  22  °C 


are  shown  in  Figs.  2  and  3,  respectively.  It  seems  that  all  candidate 
functions  work  reasonably  well  and  the  polynomial  (5)  matches  the 
OCV  points  a  little  better  for  both  cells.  The  fitting  results  at  tem¬ 
perature  =  10  °C  and  35  °C  are  similar.  It  was  argued  in  Ref.  [25]  that 
rich  features  in  low  or  high  SOC  can  be  better  fitted  by  a  high-order 
smooth  polynomial.  Therefore,  function  (5)  in  Table  2  is  selected. 

2.2.  Optimization  process  for  the  RC  circuit  and  one-state  hysteresis 
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In  Ref.  [24],  the  first-order  RC  model  with  one-state  hysteresis 
was  found  to  fit  the  LiFeP04  cell  results  the  best  while  the  first- 
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Fig.  7.  Comparison  results  of  polynomials  with  different  degrees  for  the  LiFeP04  cell. 


Fig.  8.  Comparison  results  of  the  three  types  of  analytic  functions  for  the  LiFeP04  cell. 
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are  summarized  in  Table  3.  The  same  basis  functions  are  applied  to 
the  hysteresis  parameters  as  well,  to  reduce  the  complexity  in 
optimization.  Since  each  model  parameter  is  expressed  as  an 
explicit  function  with  unknown  coefficients,  all  are  to  be  found  by 
optimization  routines.  Furthermore,  unlike  the  model  parameters 
(e.g.,  internal  resistance),  we  often  have  no  clear  knowledge  of  the 
bounds  of  the  unknown  coefficients.  The  global  multi-swarm 
particle  swarm  optimization  (MPSO)  method  [24]  is  used  to  opti¬ 
mize  these  coefficients.  However,  its  computation  load  is  very 
heavy.  If  gradient-based  methods  such  as  sequential  quadratic 
programming  (SQP)  are  used,  they  frequently  got  trapped  at  local 
minima.  Therefore,  a  two-step  optimization  procedure  is  adopted. 
The  schematic  diagram  of  the  two-step  optimization  is  shown  in 
Fig.  4.  With  loose  bounds  for  all  the  unknown  coefficients,  4000 
generations  in  the  global  MPSO  are  first  iterated  such  that  the  RMS 
error  is  less  than  30  mV.  Then,  the  3  best  solutions  (particles) 
attained  by  MPSO  are  used  as  the  initial  solutions  of  the  SQP 


algorithm.  The  optimization  is  conducted  under  each  temperature. 
Then,  we  fit  piecewise  linear  functions  to  the  optimal  solutions  for 
the  three  temperatures,  thereby  establishing  an  explicit  tempera¬ 
ture  dependency.  The  OCV  coefficients  (sixth-degree  polynomial) 
are  also  optimized  in  the  two-step  procedure,  together  with  other 
model  parameters. 

(1)  LiNMC  cell 

The  results  using  polynomial  basis  functions  of  different  degrees 
are  shown  in  Fig.  5.  The  maximum,  minimum  and  average  RMS 
errors  under  three  temperatures  are  computed.  The  hybrid  pulse 
test  dataset  is  used  as  the  training  dataset  while  the  Dynamic  Stress 
Test  (DST)  and  Federal  Urban  Dynamic  Schedule  (FUDS)  datasets 
[24]  are  used  as  validation  datasets.  The  third-degree  polynomial  is 
selected  based  on  accuracy-complexity  trade-off.  Fig.  6  shows  the 
comparison  results  of  the  third-degree  polynomial,  exponential  (2) 


Table  5 

Optimal  coefficients  of  the  analytic  functions  for  the  LiFeP04  cell. 
Parameter 
OCV 


Charging  resistance  Rq 


Discharging  Rj 


R i  (RC  network) 


Time  constant  n  (RC  network) 


Hysteresis  decaying  factor  k 


Maximum  charging  hysteresis  H 


Maximum  discharging  hysteresis  H+ 


Coefficients 

Vqc  =  0i\Z 6  +  a2 Z5  +  CH3Z4  +  d4Z3  +  OL5Z2  +  0L&  +  0i~i  (0.1  <  z  <  0.9) 
where 

ai  =  {  -0.4629(Te  -  10)  +  1.9754,  10  <  Te  <  22;  0.2506(Te  -  22)  -  3.5795,  22  <  Te  <  35} 

a2  =  {  1.0943(Te  -  10)  -  3.6208,  10  <  Te  <  22;  -0.7233(Te  -  22)  +  9.5106,  22  <  Te  <  35} 

a3  =  {  -0.8341  (Te  -  10)  +  0.0094,  10  <  Te  <  22;  0.8529(Te  -  22)  -  10.0000,  22  <  Te  <  35} 
a4  =  (0.1023(Te- 10) +  4.3001,  10  <  Te  <  22;  -0.4813(Te  -  22)  +  5.5273,  22  <  Te  <  35} 
a5  =  (0.1514(Te- 10) -3.7275,  10  <  Te  <  22;  0.1 121  (Te  -  22)  -  1.9108,  22  <  Te  <  35} 
a6  =  {  -0.0631  (Te  -  10)  +  1.3896,  10  <  Te  <  22;  -0.0049(Te  -  22)  +  0.6326,  22  <  Te  <  35} 
a7  =  {  0.0068(Te  -  10)  +  3.0629,  10  <  Te  <  22;  -0.0002(Te  -  22)  +3.1440,  22  <  Te  <  35} 
R-  =  C\Z3  +  c2z2  +  c3z  +  c4  (0.1  <  z  <  0.9) 
where 

Ci  =  {  0.0148(Te  -  10)  -  0.1670,  10  <  Te  <  22;  -0.0018(Te  -  22)  +  0.0105,  22  <  Te  <  35} 
c2  -  {  -0.0221  (Te  -  10)  +0.2458,  10  <  Te  <  22;  0.0027(Te  -  22)  -  0.0200,  22  <  Te  <  35} 
c3  -  {  0.0098(Te  -  10)  -  0.1027,  10  <  Te  <  22;  -0.0014(Te  -  22)  +  0.0149,  22  <  Te  <  35} 
c4  =  {-0.0016(Te- 10) +0.0349,  10  <  Te  <  22;  -0.0001  (Te  -  22)  +  0.0151,  22  <  Te  <  35} 
Rj  =  CiZ3  +  c2z2  +  c3z  +  c4  (0.1  <  z  <  0.9) 
where 

Cl  =  {  -0.0012(Te  -  10)  -  0.0010,  10  <  Te  <  22;  0.0106(Te  -  22)  -  0.0153,  22  <  Te  <  35} 

c2  =  (0.0005(Te- 10) +  0.0274,  10  <  Te  <  22;  -0.01 63(Te  -  22)  +  0.0333,  22  <  Te  <  35} 

c3  =  {  0.0009(Te  -  10)  -  0.0340,  10  <  Te  <  22;  0.0071  (Te  -  22)  -  0.0231,  22  <  Te  <  35} 
c4  =  {  -0.0013(Te  -  10)  +  0.0432,  10  <  Te  <  22;  -0.0014(Te  -  22)  +  0.0277,  22  <  Te  <  35} 
Ri  =  CiZ3  +  c2z2  +  c3z  +  c4  (0.1  <  z  <  0.9) 
where 

Ci  =  {  0.0272(Te  -  10)  -  0.3042,  10  <  Te  <  22;  -0.0058(Te  -  22)  +  0.0220,  22  <  Te  <  35} 
c2  =  {  -0.0424(Te  -  10)  +  0.4777,  10  <  Te  <  22;  0.0078(Te  -  22)  -  0.0309,  22  <  Te  <  35} 
c3  =  {  0.0226(Te  -  10)  -  0.2715,  10  <  Te  <  22;  -0.0022(Te  -  22)  -  0.0000,  22  <  Te  <  35} 
c4  =  {-0.0054(Te-  10) +0.0956,  10  <  Te  <  22;  -0.0007(Te  -  22)  +  0.0311,  22  <  Te  <  35} 
Ti  =  CiZ3  +  c2z2  +  c3z  +  c4  (0.1  <  z  <  0.9) 
where 

Ci  =  {0.1124(Te- 10) -4.3311,  10  <  Te  <  22;  -0.2546(Te  -  22)  -  2.9824,  22  <  Te  <  35} 
c2  =  {  -0.4379(Te  -  10)  +  1.1861,  10  <  Te  <  22;  -0.0584(Te  -  22)  -  4.0682,  22  <  Te  <  35} 
c3  =  (0.2447(Te-  10) +  6.2657,  10  <  Te  <  22;  -1.4771  (Te  -  22)  +  9.2026,  22  <  Te  <  35} 
c4  =  {  -0.9800(Te  -  10)  +  35.2236,  10  <  Te  <  22;  0.5832(Te  -  22)  +  23.4637,  22  <  Te  <  35} 
k  =  CiZ3  +  c2z2  +  c3z  +  c4  (0.1  <  z  <  0.9) 
where 

Ci  =  {  0.0379(Te  -  10)  -  0.5478,  10  <  Te  <  22;  0.0243(Te  -  22)  -  0.0925,  22  <  Te  <  35} 
c2  =  { -0.0521  (Te-  10) +0.7931,  10  <  Te  <  22;  -0.0339(Te  -  22)  +  0.1675,  22  <  Te  <  35} 
c3  =  {  0.0188(Te  -  10)  -  0.3157,  10  <  Te  <  22;  0.0124(Te  -  22)  -  0.0903,  22  <  Te  <  35} 
c4  =  {-0.0015(Te- 10) +0.0410,  10  <  Te  <  22;  -0.001 1  (Te  -  22)  +  0.0233,  22  <  Te  <  35} 
H~  =  C\Z3  +  c2z2  +  c3z  +  c4  (0.1  <z<  0.9) 
where 

Ci  =  {  -0.0248 (Te  -  10)  +  0.9525,  10  <  Te  <  22;  -0.0725(Te  -  22)  +  0.6554,  22  <  Te  <  35} 
c2  =  {  0.0306(Te  -  10)  -  1.3109,  10  <  Te  <  22;  0.1021  (Te  -  22)  -  0.9440,  22  <  Te  <  35} 
c3  =  {  -0.0085(Te  -  10)  +  0.4424,  10  <  Te  <  22;  -0.0331  (Te  -  22)  +  0.3399,  22  <  Te  <  35} 

c4  =  {  -0.0005(Te  -  10)  +  0.0060,  10  <  Te  <  22;  0.001 8(Te  -  22)  +  0.0000,  22  <  Te  <  35} 

H+  =  CiZ3  +  c2z2  +  c3z  +  c4  (0.1  <  z  <  0.9) 
where 

Ci  =  (0.0415(Te-  10) -0.0848,  10  <  Te  <  22;  -0.0091  (Te  -  22)  +  0.4134,  22  <  Te  <  35} 
c2  =  {  -0.0609(Te  -  10)  -  0.0169,  10  <  Te  <  22;  0.0165(Te  -  22)  -  0.7478,  22  <  Te  <  35} 

c3  =  {  0.0243(Te  -  10)  +  0.0944,  10  <  Te  <  22;  -0.0065(Te  -  22)  +  0.3857,  22  <  Te  <  35} 

c4  =  {  -0.0018(Te  -  10)  -  0.0392,  10  <  Te  <  22;  0.0005(Te  -  22)  -  0.0613,  22  <  Te  <  35} 
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Fig.  9.  OCV  function  of  the  LiNMC  cell. 


and  power  series  (3).  The  results  are  all  very  similar  (i.e.,  they  all 
work  reasonably  well).  The  exponential  function  shows  slightly 
better  validation  accuracy  (especially  in  FUDS  data)  and  is  quite 
simple.  Therefore,  the  exponential  function  is  used  to  describe  the 
RC  parameters  of  the  first-order  RC  model  for  the  LiNMC  cells.  The 
optimal  coefficients  of  the  analytic  functions  are  summarized  in 
Table  4. 


Table  6 

EKF  equations  for  the  SOC  estimation. 

Nonlinear  battery  model 

*k+i  =  A(zk,  Te)xk  ±  B(zk,  T e)uk  ±  o)k 

yk  =  g(xk,  4)  +  vk 

where  x,  I  and  y  are  state  variables,  current  and  voltage,  respectively. 
The  first  element  of  x  is  SOC  z.  u  is  the  model  input.  A  and  B  are  SOC 
and  temperature-dependent  matrices,  w  and  v  are  assumed  to  be 
independent,  zero-mean,  Gaussian  noise  processes 
of  covariance  matrices  Q.  and  V. 

(1)  Initialization 

Assign  the  initial  state  estimate  x0,  error  covariance  P0,  Q.and  V. 

(2)  Prediction 


I  =  A(zk_uTe)xk_,  +Biik_, 

\  Pk|k-i  =  A(zk_,Je)Pk_,A(zk_iJe)T  +  Q 

(3)  Correction 


CGk  =  f>k|k-iC(Zk|k_i,Te)T[C(zk|k_1,Te)Pk|k_1C(zk|fc_1,Te)T  +  W]  1 
\  *k  =  xk\k-i  +  GiM  ~  g(xk\k-\,Ik)] 

{Pk  =  [E  ~  GkC(zk,  Te)]Pk|k_a 
where 


C(zk|  k~\  j  Te) 


dg(*kTk), 

dXk  lxk=xm^ 


Table  7 

Matrices  of  the  state-space  models  for  the  LINMC  and  LiFeP04  cells. 
(1)  The  first-order  RC  model  of  the  LiNMC  cell 


(2)  LiFeP04  cell 

Polynomials  of  different  degrees  for  the  LiFeP04  cell  are  shown 
in  Fig.  7.  Fig.  8  shows  the  performance  of  the  third-degree  poly¬ 
nomial,  exponential,  and  power  series  models.  It  was  found  that  the 
third-degree  polynomial  function  shows  the  best  balance  between 
training/validation  accuracies  and  complexity.  The  optimized 
coefficients  of  the  analytic  functions  for  the  LiFePCU  cell  (the  first- 
order  RC  model  with  one-state  hysteresis)  are  summarized  in 
Table  5.  The  OCV  functions  of  the  two  types  of  Li-ion  batteries  are 
shown  in  Figs.  9  and  10.  It  is  obvious  that  the  OCV  of  the  LiFePCU  cell 
between  25%  and  85%  SOC  is  much  flatter  than  that  of  the  LiNMC 
cell.  The  exact  dependence  of  the  model  parameters  on  SOC  and 
temperature  are  unclear  due  to  the  complex  reactions  and 
dynamics  inside  the  cells  [26].  Furthermore,  since  we  use  lumped 
models  to  approximate  the  current— voltage  relationships,  the 
model  parameters  are  just  estimates  with  no  clear  physical 


°  exp(-±) 

-r/A  t 

*k  =  [ft  ].»(**,  Te)  = 

uk 

.B(zk,  Te)  = 

3600 Cn 

R!  - exp 

g(*kJk)  =  VOc,k-Vk-ltR0-IkR0 

where  U  is  the  voltage  across  the  RC  network  at  time  k.  At  is  the  sampling  time, 
r]  is  Coulombic  efficiency,  and  Cn  is  the  cell  capacity.  I+  and  /  are  discharging 
and  charging  currents,  respectively.  When  I  <  0  (charging),  I+  =  0,1=  /; 
when  I  >  0  (discharging),  r  =  0, 1+  =  I. 

(2)  The  first-order  RC  model  with  hysteresis  for  the  LiFeP04  cell 


zk 

*k  =  [Uk],A(zk,Te) 

K 


r1  0  0  i 

0  exp(i?)  0 

.0  0  exp(-|/c/kAt|). 


B(z/oTe) 


-Tj/S  t 
3600 Cn 
Ri  -  Riexp 
0 


0 

0 

1  -  exp(-\Klk\t\) 


SOC 


"  =  {  [h-]’7  -  0;  [L  ]  ’^  >  £(xk’k)  =  V0QM-Uk  +  hk-l+R+-I;Rj. 


meanings.  Nevertheless,  the  models  were  able  to  represent  the 
lumped  behavior  of  the  battery  well.  With  the  RMS  error  around 
merely  5  mV,  the  achieved  analytic  functions  can  be  used  to  predict 
battery  voltage  with  adequate  accuracy.  It  should  be  emphasized, 
however,  that  the  model  error  in  the  validation  dataset  is  much 
higher.  This  is  perhaps  because  the  hysteresis  model  was  trained 

Table  8 

Assumed  initial  SOC  error  of  Coulomb  counting  for  the 

LiNMC  cell. 

10  °C  22  °C  35  °C 

±3.74%  ±3.09%  ±3.28% 


Fig.  10.  OCV  function  of  the  LiFeP04  cell. 
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- EKF - Coulomb  counting  (positive  initial  error)  Coulomb  counting  (negative  initial  error) - 15% 

10 1 - 1 - 1 - 1 - 1 - 1 - 1 - 


.10 1 - 1 - 1 - 1 - 1 - 1 - 1 - 

0  0.5  1  1.5  2  25  3  35 

D  DSTtest 


C  FUOS  test 


Fig.  11.  SOC  estimation  errors  of  the  LiNMC  ceil  under  22  °C  (before  aging),  (a)  hybrid 
pulse  test;  (b)  DST  test;  (c)  FUDS  test. 


Table  9 

Assumed  initial  SOC  error  of  Coulomb  counting  for  the  LiFeP04 
cell. 


10  °c 

22  °C 

35  °C 

±6.30% 

±5.41% 

±9.65% 

using  less-transient  dynamic  dataset  and  during  the  more  signifi¬ 
cant  transient  the  model  error  increases. 

3.  EKF-based  SOC  estimator 

EKF  is  a  widely  used  state  estimation  method  for  nonlinear 
dynamical  systems.  Based  on  the  optimal  models  presented  in 
Section  2,  EKF  is  adopted  to  estimate  the  battery  SOC.  The  EKF 
equations  for  the  SOC  estimation  are  summarized  in  Table  6.  The 
matrices  for  the  state-space  models  for  the  LINMC  and  LiFeP04  cells 
are  summarized  in  Table  7. 

4.  Robustness  of  the  SOC  estimator 

To  better  evaluate  the  robustness  of  the  SOC  estimator  in  an 
environment  close  to  realistic  automotive  environment,  artificial 
Gaussian  noise  is  added  to  the  current  collected  by  the  Arbin  tester. 
The  mean  and  standard  deviation  of  the  noise  are  -0.01  A  and 
0.06  A  (around  2%  of  the  full  range),  respectively. 

4.1  Robustness  against  different  loading  profiles  and  temperatures 

The  hybrid  pulse  test,  DST  and  FUDS  datasets  collected  under 
three  different  temperatures  for  new  cells  are  used  in  this  subsec¬ 
tion.  The  data  below  20%  SOC  is  discarded,  since  the  battery  SOC  is 
seldom  permitted  to  be  less  than  20%  in  actual  driving.  To  better 


g 6 
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a>  4 


Hybrid  Pulse  DST  FUDS 

a  10°C 


|  EKF  Coulomb  counting  (positive  mital  error)  I  I  Coulomb  counting  (negative  inital  error) 


Fig.  12.  RMS  SOC  errors  under  three  temperatures  for  the  LiNMC  cell  (new  cell). 
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C  FUOS  test 

Fig.  13.  SOC  estimation  errors  of  the  LiFeP04  cell  under  22  °C  (before  aging),  (a)  hybrid 
pulse  test;  (b)  DST  test;  (c)  FUDS  test. 


evaluate  the  SOC  estimator,  its  performance  is  compared  against 
the  Coulomb  counting  method. 

(1)  LiNMC  cell 

The  OCV  is  often  used  to  determine  the  initial  SOC  of  Coulomb 
counting  using  the  OCV-SOC  relationship.  The  OCV  error  thus  leads 
to  the  initialization  error  for  Coulomb  counting.  Based  on  the  fact 
that  the  practically  available  cell-level  precision  is  around  0.1%  [31  ], 
the  OCV  error  incurred  by  voltage  sensor  is  thus  around  4  mV. 
Furthermore,  in  practice  OCV  measurement  is  subject  to  electro¬ 
magnetic  interference,  switching  of  high-current  loads  in  nearby 
wiring  and  insufficient  relaxation.  Sometimes  insufficient  relaxa¬ 
tion  can  induce  a  voltage  measurement  error  of  20  mV  [32].  Based 
on  these  considerations,  we  assume  that  the  absolute  measure¬ 
ment  error  of  the  cell  OCV  is  8  mV.  The  maximum  initial  SOC  error 
of  the  Coulomb  counting  method  can  be  algebraically  calculated 
based  on  this  assumed  OCV  error  and  the  SOC-OCV  curve.  The 
results  are  summarized  in  Table  8. 

The  initial  error  of  the  EKF  is  assumed  to  be  10%.  The  initial 
values  of  the  EKF  prediction  and  correction  gain  matrices  are  fixed 
for  all  cases.  Note  that  both  the  EKF  and  Coulomb  counting  use  the 
corrupted  current  measurement  in  their  SOC  estimation.  Fig.  11 
shows  the  SOC  estimation  results  for  22  °C.  It  can  be  seen  that 
the  EKF-based  estimator  efficiently  compensates  for  the  initial 
error  and  has  a  better  accuracy  than  Coulomb  counting  for  both  the 
training  and  validation  datasets.  The  biased  current  noise  has 
a  larger  negative  impact  to  the  Coulomb  counting  method  than  to 
the  EKF.  The  results  under  10  °C  and  35  °C  are  similar.  To  more 
clearly  compare  the  overall  performance  of  the  EKF  and  Coulomb 
counting,  the  RMS  SOC  estimation  errors  under  three  temperatures 
are  shown  in  Fig.  12.  It  is  obvious  that  the  EKF  SOC  estimator  has 
better  accuracy.  Since  the  model  has  an  inferior  accuracy  in  the 
validation  data  (DST/FUDS  dataset  excluding  the  data  below  20% 
SOC)  to  that  in  the  training  data  (hybrid  pulse  test  dataset  excluding 
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Fig.  14.  RMS  SOC  errors  under  three  temperatures  for  the  LiFeP04  cell  (new  cell). 
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Fig.  15.  Normalized  capacities  of  the  LiNMC  cell  in  the  aging  process. 


the  data  below  20%  SOC),  the  error  of  SOC  estimation  in  the  vali¬ 
dation  data  is  always  higher.  The  estimation  errors  in  the  validation 
datasets  under  three  temperatures  are  around  1%,  indicating  a  good 
robustness  against  varying  loading  profiles  and  temperatures. 

(2)  LiFeP04  cell 

The  initial  SOC  error  of  Coulomb  counting  for  the  LiFeP04  cell  is 
determined  using  the  same  method  as  the  LiNMC  cell.  Again 
assuming  a  voltage  measurement  error  of  8  mV,  the  results  for  the 
LiFeP04  cell  are  shown  in  Table  9.  The  errors  are  larger  than  those  of 


the  LiNMC  cell,  due  to  the  flatter  OCV-SOC  curves.  The  SOC  esti¬ 
mation  errors  under  22  °C  are  shown  in  Fig.  13.  The  RMS  SOC  errors 
under  three  temperatures  are  shown  in  Fig.  14.  The  results  are 
similar  to  those  of  the  LiNMC  cells.  Further,  we  can  see  that  there  is 
an  obvious  difference  in  estimation  errors  between  the  training  and 
validation  datasets.  Due  to  the  flatter  OCV,  the  SOC  estimation  of 
the  LiFeP04  cell  is  more  susceptible  to  model  error.  The  validation 
errors  are  around  2.5%,  higher  than  the  1%  accuracy  achieved  on  the 
LiNMC  cells. 

4.2.  Robustness  against  cell  aging 

The  battery  cells  always  degrade  with  time  in  real  operations. 
The  degradation  is  manifested  as  capacity  or  power  losses.  The 
capacity  loss  is  partially  due  to  loss  of  cyclable  Li  Ions  caused  by 
multiple  factors,  such  as  cathode  structure  degradation,  side  reac¬ 
tions,  passivation  form  and  lithium  plating  at  the  anode  [33].  The 
power  loss  is  mainly  due  to  the  increase  in  cell  resistance.  It  is 
important  to  study  the  robustness  of  the  SOC  estimator  since 
robustness  against  cell  aging  is  of  great  significance  for  BMS  of 
electrified  vehicles. 

(1)  LiNMC  cell 

Fig.  15  shows  the  calibrated  capacities  under  three  temperatures 
in  the  degradation  process  of  the  LiNMC  cell  at  22  °C.  Flere,  we  used 
data  up  to  cell  age  of  342  cycles  at  22  °C  to  analyze  the  robustness  of 
the  EKF  estimator.  The  capacity  reduced  by  about  10%  at  10  °C  when 
the  cell  is  342  cycles  old  at  22  °C— which  is  roughly  half  way 
through  the  typical  life  expectancy  accepted  for  automotive  appli¬ 
cations  (a  battery  pack  is  determined  to  be  “at  the  end  of  its  life” 
when  its  capacity  loss  reaches  20%).  The  three  characterization  test 
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Fig.  16.  Comparison  results  of  a  sequence  of  estimators  for  the  LiNMC  cell  (after  aging)  with  different  level  of  model  parameter  update. 
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Fig.  17.  Comparison  results  of  a  sequence  of  estimators  for  the  LiFeP04  cell  (after  aging)  with  different  level  of  model  parameter  update. 


datasets  right  after  342  cycles  are  used  to  evaluate  the  SOC  esti¬ 
mator.  The  initial  SOC  error  and  tuning  parameters  of  the  EKF  are 
kept  unchanged  when  we  compare  a  sequence  of  estimators  with 
increasing  level  of  parameter  updates.  Fig.  16  shows  the  compar¬ 
ison  results.  It  can  be  seen  that  if  no  parameters,  including  the 
capacity  Cn,  are  updated,  the  estimator  performance  deteriorates 
significantly.  The  RMS  errors  of  the  SOC  estimation  are  clearly 
larger  than  those  shown  in  Fig.  12.  Updating  internal  ohmic  resis¬ 
tance  Ro  can  improve  the  SOC  estimation  in  most  cases.  Updating 
parameters  of  the  RC  network  in  addition  to  Ro  does  not  help  much. 
If  we  update  all  model  parameters,  including  the  cell  capacity  Cn, 
the  RMS  error  can  be  reduced  to  below  1%  in  most  cases. 

It  should  be  noted  that  updating  some  of  the  model  parameters 
is  quite  easy  to  do.  For  example,  the  internal  ohmic  resistance  Ro 
can  be  easily  obtained  even  under  simple  pulse  loads.  This  simple 
update  reduces  model  prediction  error  by  more  than  half.  Updating 
RC  parameters  or  the  cell  capacity  is  harder  onboard.  It  may  be 
possible  to  obtain  those  model  parameters  during  charging  process, 
when  the  charging  current  can  be  manipulated  and  more  compu¬ 
tation  power  is  available. 

(2)  LiFeP04  cell 

The  LiFeP04  cells  we  tested  show  much  slower  degradation.  We 
have  finished  2600  cycles  and  the  capacity  loss  at  22  °C  is  only 
around  3%.  The  three  characterization  test  datasets  collected  right 
after  2587  cycles  are  used.  The  RMS  errors  of  the  estimators  with 
different  levels  of  model  parameter  updates  are  shown  in  Fig.  17. 
Again,  the  estimator  without  parameter  update  shows  poor  per¬ 
formance— this  is  because  while  the  cell  capacity  does  not  reduce 


significantly,  its  dynamic  behavior  changes  noticeably.  The  RMS 
errors  in  the  training  datasets  (hybrid  pulse  test)  are  worse  than  the 
LiNMC  cell,  due  to  the  very  flat  OCV-SOC  curve  of  the  LiFeP04  cell. 
The  RMS  errors  in  the  validation  datasets  (DST  and  FUDS)  are 
similar  to  the  LiNMC  cell.  The  nominal  model  of  the  LiFeP04  cell  is 
less  accurate  after  degradation.  Also,  updating  internal  resistance 
Ro  improves  the  robustness  of  the  SOC  estimator  noticeably. 
Updating  the  RC  network  and  capacity  does  not  help  much,  since 
capacity  loss  is  small.  If  all  the  model  parameters  are  updated,  the 
SOC  estimator  has  the  best  robustness  based  on  the  overall  evalu¬ 
ation  of  all  the  9  datasets.  In  practice,  however,  it  is  challenging  to 
update  all  the  model  parameters.  Considering  both  the  model 
accuracy  and  requirement  for  parameter  update,  we  concluded  that 
updating  the  internal  resistance  is  the  most  cost-effective  update 
one  can  use  to  improve  the  robustness  of  SOC  estimation  method 
against  cell  aging. 

5.  Conclusions 

The  robustness  of  SOC  estimation  algorithms  of  two  types  of  Li- 
ion  batteries  is  studied  against  varying  loading  profiles,  tempera¬ 
tures  and  aging  levels.  The  preferred  model  structures  determined 
in  our  previous  work  [24]  for  the  LiNMC  and  LiFeP04  cells  are  used. 
The  optimal  functions  of  SOC  and  temperature  explicitly  describing 
the  model  parameters  are  selected  from  commonly  used  generic 
functions  for  battery  modeling.  The  sixth-degree  polynomial  is 
found  to  depict  OCVs  of  both  cells  well.  The  exponential  function  is 
determined  to  be  the  best  for  the  RC  parameters  of  the  first-order 
RC  model  for  the  LiNMC  cell;  the  third-degree  polynomial  is 
preferred  for  the  RC  and  hysteresis  parameters  of  the  first-order  RC 
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model  with  hysteresis  for  the  LiFeP04  cell.  The  coefficients  of  these 
generic  functions  are  optimized  by  a  combined  MPSO  and  SQP 
procedure.  Based  on  the  optimal  model  parameters,  the  Extended 
Kalman  Filter  technique  is  used  to  estimate  the  battery  SOC.  Given 
the  training  datasets  from  new  cells,  the  robustness  of  the  SOC 
estimator  against  varying  loading  profiles  and  temperatures  is 
evaluated  and  compared  against  the  Coulomb  counting  method. 
The  EKF-based  estimator  performs  better  than  the  Coulomb 
counting  method,  and  the  RMS  estimation  errors  under  validation 
datasets  under  three  temperatures  are  around  1%  and  2.5%  for  the 
LiNMC  and  LiFeP04  cells,  respectively.  Due  to  a  flatter  OCV-SOC 
surface,  the  SOC  estimation  of  the  LiFeP04  cell  is  more  vulnerable  to 
cell  aging  and  if  the  model  is  not  updated,  its  performance  was 
found  to  be  unacceptable.  Updating  the  internal  ohmic  resistance  of 
the  nominal  model  improves  the  performance  of  the  estimator 
significantly  (by  about  half).  Updating  all  the  model  parameters 
achieves  even  better  performance  however  this  may  not  be  easy  to 
achieve  in  real-time  driving. 
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